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CHAPTER  I 

INTRODUCTION 

New  system  requirements  are  more  and  more  being  directed  toward  the 
attainment  of  adequate  overall  system  characteristics  such  as  reduction 
of  maintenance  costs  and  general  support  as  well  as  high  performance 
requirements.  As  system  requirements  become  more  Involved  and  more 
complex,  it  becomes  increasingly  necessary  to  speak  of  system  charac¬ 
teristics  in  definite  quantitative  terms. 

In  recent  years  such  terms  as  capability,  dependability  and  cost 
effectiveness  have  been  quantified  In  an  attempt  to  describe  system 
performance.  However,  the  most  common  term  used,  and  the  one  most 
specified  in  government  contracts  is  inherent  availability, 

Blanchard  and  Lowery  [2]  define  inherent  availability  as  "  .  .  , 
the  probability  that  a  system  or  equipment,  when  used  under  stated  condi¬ 
tions,  without  consideration  for  any  scheduled  or  preventive  cction,  in 
an  ideal  support  environment  (ie,,  available  tools,  spares,  manpower, 
data,  etc.),  shall  operate  satisfactorily  at  a  given  point  in  time.  It 
excludes  ready  time,  preventive  maintenance  downtime,  logistics  time, 
and  waiting  or  administrative  downtime,'*  Inherent  availability  can  be 
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paper. 
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expressed  as 

.  mtbf 

AI  MTBF  +  MTTR 

Hare  MTBF  la  the  mean  time  between  failures  and  MTTR  le  the  mean  time 
to  repair. 

In  order  to  co-ordinate  the  efforts  of  different  groups  con¬ 
cerned  with  different  system  characteristics,  and  to  eliminate  the 
hazards  of  guesswork  in  achieving  the  overall  system  requirement,  it 
is  necessary  to  establish  a  procedure  for  determining  the  detailed 
specifications  for  the  various  components  that  make  up  a  system.  The 
process  of  assigning  availability  parameters  to  individual  components 
to  insure  the  attainment  of  the  system  availability  goal  is  termed 
availability  allocation. 

The  problem  that  is  solved  in  thiB  paper  assumes  that  an  avail¬ 
ability  requirement  has  been  determined  based  on  operational  need,  and 
the  paper  is  concerned  with  allocation  of  parameters  to  achieve  that 
requirement . 

The  allocation  of  system  availability  involves  solving  the  inequal¬ 
ity 

f(M,  ^2,  ...  An)  >  A* 

where  A*  is  the  system  availability  requirement,  Ai  is  the  allocated 
availability  for  subsystem  i,  and  f  is  the  functional  relationship 
between  unit  and  system  availability.*  The  above  inequality  has  an 
infinite  number  of  solutions  if  no  restrictions  are  placed  on  the  allo¬ 
cation.  The  problem  is  to  establish  a  solution  procedure  which  yields 

• 

All  availabilities  referred  to  throughout  the  paper  are  assumed  to  be 
Inherent  availabilities. 
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a  unique  or  limited  number  of  solutions.  The  approach  used  in  this 
paper  is  to  minimize  the  system  life  cycle  cost  subject  to  the  con¬ 
straint.  of  meeting  the  system  availability  goal  by  the  method  of 
Lagrange  multipliers.  Mean  repair  times  and  failure  rates  are  first 
allocated  on  a  minimum  cost  basis  with  the  constraint  inactive,  and 
second,  with  the  constraint  active. 

The  solution  of  this  problem  for  the  availability  parameters  MTBF 
and  MTTR  is  accomplished  under  the  assumptions  of  constant  failure  rates 
(exponential  time-to-failure  distributions)  and  independently  failing 
components.  Also,  a  series  configuration  is  assumed  in  which  all  sub¬ 
systems  are  necessary  for  the  system  to  function.  The  limitations  of 
these  assumptions  are  discussed  in  Chapter  III, 

A  general  solution  procedure  is  followed  to  +he  extent  feasible. 
Therefore,  a  solution  set  of  equations  is  expressed  under  the  assump¬ 
tions  listed  above.  This  set  of  equations  may  be  used  with  any  appli¬ 
cable  cost  functions.  Specific  cost  functions  are  considered  for 
detailed  analysis  and  solution  in  this  paper.  However,  by  substitu¬ 
ting  any  differentiable  cost  functions  into  the  solution  equations  of 
Chapter  III,  a  set  of  2n  +  1  equations  in  2n  +  1  unknowns  could  be 
obtained,  where  n  is  the  number  of  subsystems  being  considered.  One 
is  still  faced  with  the  problem  of  solving  the  set  of  simultaneous 
equations  for  the  allocated  failure  rates  and  repair  times,  but  this 
presents  no  serious  drawback  since  computerized  numerical  techniques 
are  available  for  the  solution  of  such  systems  of  equations.  The 
methods  used  in  this  paper  should  be  helpful  in  providing  solution 
techniques  for  similar  cases, 

^  _ _ _ _ _ _ 


The  costs  considered  in  this  paper  axe  research  and  development, 
production  and  procurement,  sparing,  and  maintenance  costs.  Those 
elements  of  life  oyole  cost  not  dependent  on  failure  rate  or  repair 
time  are  not  considered  in  this  paper,  while  it  is  assumed  that  other 
costs  not  mentioned  above  which  may  be  dependent  on  failure  rate  or 
repair  time  can  be  included  in  one  of  the  above  categories.  This  will 
be  discussed  further  in  Chapter  IV. 

Development  and  production  costs  are  assumed  negatively  correlated 
with  both  failure  rate  and  repair  time.  That  is,  high  costs  are  asso¬ 
ciated  with  low  failure  rates  or  repair  times  while  relatively  low  costs 
are  associated  with  high  failure  rates  and  repair  times.  The  reason  for 
separating  development  cost  and  production  cost  will  be  discussed  in 
Chapter  III,  Sparing  costs  are  assumed  to  be  proportional  to  the  pro¬ 
curement  cost  of  spares  times  the  expected  number  of  failures  over  the 
life  cycle  whil«»  the  maintenance  cost  is  assumed  to  have  one  component 
proportional  to  the  number  of  hours  of  maintenance  required  over  the 
life  cycle,  and  another  component  proportional  to  the  number  of  failures 
which  occur  during  the  life  cycle.  The  development  of  this  model  is 
given  in  Chapter  IV. 

There  are  basically  two  situations  in  which  it  would  be  desirable 
to  allocate  availability.  The  first  case  is  that  during  the  conceptual 
stage  of  system  development  when  no  design  has  been  formalized,  an 
allocation  based  on  estimated  costs  of  reliability  and  maintainability 
development  can  give  the  optimal  set  of  component  parameters  that 
should  be  established  aa  system  goals.  Second,  is  the  case  where  a 
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particular  design  has  been  proposed  or  developed  that  falls  short  of 
the  system  availability  requirement.  The  solution  given  in  this  paper 
is  applicable  to  either  case. 

In  Chapter  III  the  general  allocation  technique  is  described  and 
limitations  are  discussed.  The  specific  coot  functions  are  developed 
in  Chapter  IV  and  the  solution  to  this  case  is  given  in  the  form  of  a 
computer  program  in  the  Appendix.  In  Chapter  V  an  example  problem  is 
stated  and  solved  and  application  of  the  computer  program  is  discussed. 
In  Chapter  VI  some  estimates  of  the  sensitivity  of  failure  rates  and 
repair  times  to  the  various  cost  factors  are  obtained  by  repeated  com¬ 
puter  runs.  The  summary,  conclusions,  and  recommendations  for 
further  study  are  given  in  Chapter  VII, 

A  brief  literature  survey  is  given  in  Chapter  II,  discussing 
pertinent  references  to  reliability-maintainability  allocation  and  life 
cycle  cost  models. 


CHAPTER  II 


LITERATURE  REVIEW 

Although  a  number  of  reliability  allocation  techniques  exist, 
rather  limited  work  has  been  carried  out  in  the  field  of  maintainability 
allocation  and  reliability-maintainability  allocation.  The  obvious 
reason  for  this  discrepancy  is  that  there  simply  has  not  been  enough 
time  in  the  brief  history  of  maintainability  for  these  developments  to 
occur. 

A  brief  description  of  the  current  reliability-maintainability 
trade-off  and  allocation  techniques  will  help  establish  a  background 
for  the  development  of  the  paper. 

A  basic  allocation  technique  is  given  in  the  Martin-Marietta  pre¬ 
pared  AMCP-?05-t  [6]  which  calculates  repair  time  by  summing  each  com¬ 
ponent  repair  time  multiplied  by  its  expected  fraction  of  occurrence, 

MTTRS  -  Jt  [MTTRi|-^-] 

1=1 

where  MTTRS  is  the  mean  system  repair  time,  MTTR^  is  the  mean  repair 
time  for  subsystem  i,  and  is  the  failure  rate  of  subsystem  i.  Thus, 
given  a  system  availability  requirement  and  subsystem  failure  rates,  the 
system  repair  time  goal  can  be  determined.  If  estimated  subsystem 
repair  times  are  inadequate  to  meet  the  availability  requirement,  then 
general  recommendations  are  given  as  to  which  subsystems  or  components 
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to  improve.  Ideally,  subsystems  exhibiting  the  higher  failure  rates 
require  a  low  repair  time,  while  those  with  low  failure  rates  can  have 
higher  repair  times.  However,  no  quantitative  relationships  have  been 
developed,  and  no  attempt  is  made  to  determine  the  effect  of  the  changes 
on  system  ccs*  _ 

Eimstaa  £3]  presents  an  availability  allocation  method  which  is 
based  on  optimizing  the  cost  effectiveness  ratio  with  checks  So  that 
repair  time  and  failure  rate  extremes  are  not  exceeded.  He  also  pre¬ 
sents  a  design  trade-off  analysis  which  is  graphically  solved  for 
specific  effectiveness  models. 

Blanchard  and  Lowery  [2]  present  a  reliability-maintainability 
trade-off  technique  which  consists  of  evaluating  the  total  system  cost 
savings  for  each  of  a  limited  number  of  design  alternatives. 

Goldman  and  Slattery  present  a  reliability-maintainability 
allocation  method  )n  which  isoavailability  functions  are  defined  and 
failure  rates  and  repair  times  are  allocated  on  the  basis  of  incre¬ 
mental  cost  analysis.  However,  no  provision  is  made  for  allocating  a 
specific  availability  to  a  subsystem, 

Nye  [8j  has  developed  a  cost  based  allocation  technique  which 
considers  maintenance,  sparing,  production  and  development  costs}  how¬ 
ever,  it  is  not  possible  to  allocate  failure  rates  and  repair  times 
simultaneously.  Either  failure  rates  are  fixed  and  repair  times  allo¬ 
cated  by  a  trial  and  error  solution  procedure  or  repair  times  are  fixed 
and  failure  rates  are  allocated. 


Messer  [7]  has  attacked  the  problem  by  assuming  that  differentiable 
cost  functions  exist  or  can  be  defined  which  give  the  changes  in  system 


8 


development  costs  as  functions  of  repair  times  and  failure  rates.  He 
then  minimizes  this  cost  function  subject  to  an  availability  constraint. 
This  paper  is  basically  an  extension  of  Messer's  paper  to  Include  pro¬ 
duction,  sparing  and  maintenance  costs. 

Although  numerous  articles  have  been  written  on  life  cycle  costing 
in  recent  years,  relatively  few  life  cycle  cost  modexs  have  been  pro¬ 
posed.  A  recent  paper  by  Koch  £5]  defines  basic  cost  categories  that 
must  be  included  in  such  a  model  and  develops  a  model  from  which  the  life 
cycle  cost  of  fixed  place  equipment  can  be  obtained.  However,  as 
stated  before,  this  paper  is  concerned  only  with  those  elements  of  life 
cycle  cost  which  are  dependent  on  repair  times  and  failure  rates. 

A  general  model  for  the  cost  based  allocation  problem  is  developed 
in  Chapter  III  in  which  the  method  of  Lagrange  multipliers  is  used  to 
formulate  the  solution. 
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CHAPTER  III 

DEVELOPMENT  OF  THE  MODEL 


As  previously  stated  the  inherent  availability  of  any  system  can 
be  expressed  as 


„  mtbf 

AI  0  MTBF  +  MTTR. 


(1) 


If  one  considers  a  series  system  made  up  of  components  having  constant 
failure  rates,  then  the  system  failure  rate  can  be  written  as  the  sum 
of  the  individual  failure  rates  if  all  components  fail  independently 
[l].  Thus, 

XEeA}+A.2  +  ,,,+  ^n 
and 

MTBF  "  ~jr -  (3) 

iMi 

where  As  is  the  system  failure  rate  and  Ap  is  the  failure  rate  of  sub¬ 
system  i. 

The  system  repaii  time  can  be  thought  of  as  the  weighted  mean  of 
the  components*  repair  times  (Mp)  where  the  weighting  factors  are  the 
ratios  of  the  component  failure  rates  to  the  total  failure  rate.  Thus, 
if  each  component’s  repair  time  is  multiplied  by  its  expected  fraction 
of  occurrence  and  the  results  are  summed  the  expected  system  repair  time 
is  obtained  as 
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Security  Classification 


Xj  X>  Xj) 

MTTR  -  ■j4ll  +  ^2  +  ... 

1  n  1  n 

“  1^1  ■  ~Y~T[  i-l 

1-1 

It  should  be  noted  here  that  the  expression  for  MTTR  assumes  that  all 
M< 's  are  mutually  exclusive  which  implies  that  when  a  failure  occurs  in 
one  subsystem  all  subsystems  are  turned  off  so  that  no  other  failure 
can  occur  while  the  system  is  being  repaired. 

Substituting  the  expressions  for  MTBF  and  MTTR  into  Equation  1 
yields 

A-  1  (5) 

A  1+2  XjMi 

for  a  system  in  which  all  subsystems  are  necessary  for  proper  system 
operation,  failures  occur  Independently  of  each  other  and  all  compo¬ 
nents  exhibit  constant  failure  rates.  The  limitations  of  these  assump¬ 
tions  are  discussed  below. 

The  assumption  of  constant  failure  rates  or  equivalently,  an  expo¬ 
nential  time  to  failure  density,  dictates  that  the  failure  rate  be  con¬ 
stant  throughout  the  life  of  the  equipment.  In  addition  to  being 
necessary  to  solve  analytically  many  system  relationships,  most  elec¬ 
tronic  equipment  appears  to  fit  this  assumption  reasonably  well  for  a 
long  period  of  time  after  the  burn  in  stage.  For  systems  that  exhibit 
weao.'-out,  or  increasing  failure  rate,  well  designed  and  implemented 
preventive  maintenance  policies  will  reduce  or  eliminate  systematic 
failures  and  provide  systems  that  exhibit  fairly  constant  failure  rates 

*A11  summations  are  from  1  to  n  unless  otherwise  stated. 
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over  long  periods  cf  time.  Thus,  It  seems  reasonable  that  the  expo¬ 
nential  time  to  failure  density  would  not  only  be  a  good,  relatively 
straightforward  approximation  for  a  system,  but  would  probably  be  quite 
realistic  when  considering  the  total  life  of  a  large  group  of  operating 
systems . 

The  series  assumption  is  somewhat  restrictive  in  that  no  provision 
can  v:  made  for  operation  in  a  degraded  mode.  By  judicious  definition 
of  components  or  units  that  make  up  a  subsystem  one  could  obtain 
reasonable  approximations  to  a  series  sys+^m,  ■'"or  example,  two  units  in 
parallel  could  be  combined  as  one  unit  using  the  combined  failuie  rates 
as  the  resultant  unit  failure  rate.  It  should  be  pointed  out,  however, 
that  a  parallel  combination  of  elements  with  constant  failure  rates 
does  not  result  in  a  unit  with  a  constant  failure  rate,  and  therefore, 
the  constant  failure  rate  assumption  must  be  reanalyzed. 

The  assumption  of  independently  failing  subsystems  or  units  is 
usually  mede  in  reliability  discussions ,  This  assumption  is  subject 
to  some  question  since,  in  most  mechanical  or  electrical  systems, 
failure  of  one  component  can  damage  or  degrade  other  components.  How¬ 
ever,  by  thorough  failure  mode  analysis  some  of  the  difficulty  could  be 
removed  by  combining  in  the  mean  time  to  reDair  the  time  required  to 
repair  all  items  which  failed  or  were  damaged  by  the  initial  chance 
failure.  This  assumption,  like  the  others,  however,  does  somewhat 
restrict  application  of  the  model  developed  in  this  paper.  It  is  also 
assumed  that  when  a  failure  occurs  the  entire  system  is  shut  off  while 
the  failed  subsystem  is  being  repaired.  The  reasons  for  this  restriction 
are  that  It  simplifies  formulation  of  the  availability  constraint,  and 
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it  enable.'  determination  of  the  total  operating  time  for  each  &uu- 
system.  If  each  subsystem  operates  the  same  amount  of  time,  determina¬ 
tion  of  the  total  operating  time  is  straightforward  when  there  is  an 
availability  constraint.  However,  if  this  assumption  is  not  made  the 
operating  time  for  subsystem  i  is  a  function  of  Xy  and  The  effect 

of  ignoring  this  assumption  would  also  result  in  an  erroneoi  .  expres¬ 
sion  for  availability,  Equation  5*  since  it  would  be  equivalent  to 
ignoring  the  probability  that  two  or  more  subsystems  are  being  repaired 
at  the  same  time.  Therefore,  this  assumption  also  rest  lets  applica¬ 
bility  of  the  model  developed  In  this  paper,  and  full  consideration 
should  be  given  to  the  assumptions  whenever  the  model  is  applied  to 
an  allocation  problem. 

The  allocation  of  availability  parameters  to  the  various  units  or 
subsystems  of  a  complex  system  involves  making  decisions  as  to  which 
components  to  consider  for  availability  improvement  and  to  w-lch 
characteristics,  repair  times  and/or  failure  rates,  improvement  effort 
should  be  r.pplied.  The  remainder  of  this  chapter  discusses  the  basic 
development  of  these  decision  criteria. 

For  cases  l  and  2  of  the  following  development  it  is  assumed  that 
a  design  has  been  proposed  that  falls  short  of  the  system  availability 
goal.  That  is  A  <  A',  where  A'  Is  the  system  availability  requirement. 
In  case  1  it  is  assumed  that  the  reliability  of  the  system  is  at  its 
highest  achievable  level  and  only  repair  times  can  be  improved.  Then, 
in  case  2  the  opposite  is  considered,  that  is,  only  the  failure  rates 
can  be  improved.  In  case  3  the  first  2  cases  are  combined  to  give  the 
overall  optimal  set  of  improvements. 
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The  components  of  life  cycle  c  st  considered  in  this  paper  are 

1)  df^(XjL),  where  df-p  defines  the  functional  relationship 
for  the  research  and  development  costs  of  changing  the 
failure  rate  Xi,* 

2)  Pf  i ( Xj_ )  a  where  pp^  defines  the  functional  relationship 
for  the  production  and  procurement  costs  of  changing  the 
failure  rate  X^f 

3)  dmp(Mp),  which  represents  the  research  and  development 
costs  associated  with  changing  repair  time  Mi* 

4)  which  represents  the  production  and  procurement 
costs  associated  with  changing  repair  time  Mi, 

5)  si(X1,  Mi),  the  sparing  costs  associated  with  subsystems 
i,  and 

6)  f^(X^,  Mi),  +he  maintenance  costs  associated  with 
subsystem  i, 

"Here,  the  cost  of  maintenance  manuals  and  test  equipment  should  be 
included  in  production  costs  if  they  are  dependent  on  repair  times  and 
failure  rates.  Other  comprnents  of  li^e  cycle  cost  exist  which  may  be 
dependent  on  repair  times  and  failure  rates;  however,  it  is  felt  that 
they  can  be  included  in  the  above  costs  with  no  loss  of  generality  to 
the  model.  However,  one  question  arises  cis  to  whether  the  operating 
costs  should  be  Included.  Operating  costs,  per  se,  are  relatively 
independent  of  repair  times  and  failure  rates.  The  cost  of  "lost 
operation",  or  decline  in  profits  due  to  downtime,  can  be  included  in 

* 

It  is  assumed  here  that  it  is  possible  to  change  the  failure  rate  of  a 
subsystem  without  affecting  its  repair  time  and  vice  versa. 
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maintenance  costs,  but  it  should  not  be  included  if  the  total  system 
operating  time  is  fixed.  In  other  words,  if  a  system  is  designed  for 
a  ten  year  life  span  with  a  .95  availability  the  total  operating  time 
is  made  independent  of  the  individual  repair  times  and  failure  rates. 
Case  ii  Reliability  Fixed  at  the  Highest  Achievable  Level 
From  Equation  5 


i  ^  At 

where  and  represent  the  presently  achieved  values  of 
rate  i  and  repair  time  i  respectively. 

Let  be  the  allocated  repair  time  for  subsystem  i  so  that 


(6) 

failure 


- i - A', 

i  +  Z*iMi  (7) 

Rewriting,  Equation  7  yields 

LAiMi  “  ~  !•  (8) 

Using  the  cost  components  previously  defined  the  total  cost*  becomes 

T.C,  -  +  Pfl  (*i)  +  *mi  (Mi)  +  Pmi  0 

+  Si(^itMi)  +  fi(^i,Mi)].  (9) 

No  definite  functional  relationships  for  the  cost  components  are  con¬ 
sidered  in  this  chapter. 

Thus,  the  allocation  problem  is  expressed  as  the  problem  of  mini¬ 
mizing  the  total  cost  of  availability  improvement  (in  this  case,  repair 
time  improvement)  subject  to  the  constraint  of  meeting  the  system 


This  does  not  represent  the  total  life  cycle  cost,  but  that  portion  of 
It  which  is  of  interest  in  allocating  repair  times  and  failure  rates. 


availability  requirement.  Thus,  the  problem  can  be  stated  as 
minimize  T.C.  -  Ejd^C^)  +  PfiC^)  +  <1*1(1%)  +  p*i(Mi) 

+  Bi(%*Mi)  +  fi(%*Mi)] 

subject  to 

^iMi  mAr~U 

If  the  constraint  equation  is  written  to  be  equal  to  zero,  that  is 

£%Mi  -  pr  +  1-0,  (10) 

the  Lagrangian  can  be  written 

A  -  E[dfi(^!)  +  Pfl(ii)  +  daiCMi)  +  PmiCMi)  -i  8i(£i#Mi) 

+  fi(^i,Mi)]  +  Q[£(tlMi)  -  +  l]  (11) 

where  Q  is  the  Lagrange  multiplier.  Taking  derivatives  with  respect  to 
Mi  and  Q  and  setting  them  equal  to  zero  yields 

<)  A  „  <Jdml(Hl)  jPmj(Hi)  aBj(^i»«i) 
dMi  a  Mi  aMi  a  mi 

+  iMkiM  ♦  <&.  0  (lz) 

a  Mi 

for  i  ■  1,  2,  ,.,  n,  and 

3A  lXiMi  -  Jr  +  1  -  0,  (13) 

dQ 

The  above  set  of  equations  contains  n  +  1  unknowns  Ml,  M2,  *,,, 

Mq,  and  Q.  Since  there  are  n  +  1  equations  an  unique  solution  could  be 
obtained  if  definite  forms  of  the  cost  functions  are  inserted. 

Case  2i  Maintainability  Fixed  at  the  Highest  Achievable  Level 

In  this  case  the  availability  requirement  is  expressed  as  a  func¬ 


tion  of  the  achieved  repair  times  and  the  allocated  failure  rates. 
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1  +  LXifli  "  A  (i4) 

where  X^  is  the  allocated  failure  rate  for  subsystem  i.  The  Lagrangian 
can  now  be  written  as 

A  ”  +  Pfi(^-i)  +si(^li^i)  +  fi(^i*l?i)l 

+  Q(£X^I?i  -  "jr  ^)*  (15) 


Differentiating  with  respect  to  X^  and  Q  yields 


<) A  „  ^dfi(Xj)  Jpfl(Xj)  ^8j(Xi  A) 

j  Xi  j  xi  5x1  rxi  “ 


T^i 


+  A 


0 


(16) 


for  L  -  1,  2,  ...»  n,  and 


^4  -  LXi!?i  -  i_  +  1  -  0.  (1?) 

Jk*  A 

The  above  set  of  equations  also  contains  n  +1  unknowns  Xj ,  X2, 

Xn  and  Q,  and  therefore  could  be  solved  if  definite  forms  of  the 
functions  are  inserted.  The  third  case  is  considered  next  in  wnich 
both  failure  rates  and  repair  times  can  be  improved. 

Case  3i  Repair  Times  and  Failure  Rates  can  be  Improved 

The  availability  requirement,  expressed  as  a  function  of  both 
allocated  repair  times  and  allocated  failure  rates,  becomes 


A' 


1 

1  +  LX^* 


(18) 


The  total  cost  is  now  expressed  ast 


T.C.  “  lOifiCXi)  +  Pfi(Xt)  +  dml(Mi)  +  pmi(Mi)  +  SiCXi.Mi)  +  fiCXi.Mi)] 


and  the  Lagrangian  for  the  general  case  becomes 

^  “  ^dfi(*i)  +  Pfi(^l)  +  dmi(Mi)  +  p,ni(Mi)  +  s^X^Mi) 

+  QClXXiMi)  -  p-  +l]. 

Differentiating  with  respect  to  X^,  and  Q  yields 


AA  -  ^dfi(Xl)  +  js1(Xi,Mi) 

jXi  "<)  Xi  jTi  Dd 


^fi(Xi.Mi) 

+  - pp—  +  QMi  -  0. 


for  i  -  1,  2,  n 


^ A  .  Jdmi(Mi)  ^  ^Pmi(Mi)  +  Js^Xi.Mi) 
JhT  )  Ml  JUT”  4  Mi 

+  ~  +  Q^i  “  o. 

d 


for  1  c  If  2 9  •••«  n 


jq  -  EX^  -jt+1-0,  (23) 

The  above  set  of  equations  contains  2n  +  1  unknowns  Xl ,  X2,  ,,,, 
Xn»  Mx,  H2i  mu  Mn  and  Q,  Since  there  are  2n  +  1  equations,  a  unique 
solution  should  be  obtained  if  definite  forms  of  the  cost  functions  are 
inserted.  It  should  be  pointed  out  that  this  set  of  equations  is 
general  in  that  any  cet  of  differentiable  cost  functions  could  be 
Inserted  and  solutions  obtained.  However,  one  should  keep  in  m'nd  the 
assumptions  involved  in  the  development. 


In  Chapter  IV  specific  cost  functions  are  discussed  and  analyzed. 
The  method  of  solution  used  may  provide  guide  lines  for  solutions  to 
other  cost  functions. 


CHAPTER  IV 


A  SOLUTION  TO  THE  ALLOCATION  PROBLEM 

In  this  chapter  cost  considerations  associated  with  availability 
allocation  are  discussed,  and  specific  cost  functions  are  considered 
for  solution  according  to  the  set  of  solution  equations  obtained  in 
Chapter  III. 

The  objective  of  an  allocation  program  should  be  to  allocate 
system  parameters  to  the  individual  units  or  subsystems  In  such  a  way 
as  to  minimize  the  costs  of  owning  a  system  while  meeting  system  opera¬ 
tional  performance  requirements.  For  an  availability  allocation  problem, 
the  costs  of  concern  are  those  associated  with  failure  rates  and 
repair  times.  These  cos+s  can  generally  be  assigned  to  the  following 
categories i 

1.  Development  costs, 

2.  Production  costs, 

3.  Inventory  costs, 

4.  Downtime  costs. 

As  stated  in  Chapter  I,  this  paper  is  not  concerned  with  establish¬ 
ing  an  availability  requirement.  It  is  assumed  that  an  availability  goal 
has  been  established  and  the  problem  is  to  achieve  that  goal  at  minimum 
cost.  Function'  describing  the  production  costs  of  a  system  versus  fail¬ 
ure  rate  or  repair  time  should  have  the  characteristics  of  the  curve 
shown  in  Figure  1, 
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Production  Cost 


~  "  '  or 

A^  allocated  failure  rate 
Mi  allocated  repair  time 

Figure  i.  production  Cost  Versus  Repair  Time  or  Failure  Rate 


A  function  of  this  type  suggests  that  production  costs  vary 
inversely  with  failure  rate  and  repair  time  or  that  high  cost3  are 
associated  with  low  failure  rates  and  repair  times  and  low  costs  with 
high  failure  rates  and  repair  times,  and  that  the  cost  of  producing  a 
unit  with  zero  repair  time  or  failure  rate  is  infinite  Bince  it  is 
impossible. 

The  function  used  to  describe  production  costs  for  this  paper  is 


where 


Pi(Ai#Mi) 


(py)i  Wi 

Al  Mi 


(23) 


(1)  Pi(^ifMi)  1®  the  cost  in  dollars  of  producing  subsystem  1  with 
repair  time  Mi  and  failure  rate  Ai, 

(2)  (PF)i  is  the  cost  factor  associated  with  the  difficulty  of 
producing  subsystem  i  with  failure  rate  Ai,  and 

(3)  (FM)i  is  the  cost  factor  associated  with  the  difficulty  of  pro¬ 
ducing  subsystem  i  with  repair  time  Mi* 

The  cost  associated  with  downtime  can  be  separated  into  two  general 


categories.  The  first  is  that  which  is  proportional  to  the  length  of 


tine  the  system  is  down.  This  term  would  include  such  factors  as  lost 
output  during  repair  time  and  manpower  costs  to  effect  repairs.  The 
second  term  can  be  described  as  a  maintenance  set  up  oost  whioh  is 
proportional  to  total  number  of  system  failures,  which  would  include 
system  shut  down  and  start  up  costs  and  other  factors  not  dependent  on 
the  duration  of  system  downtime.  Thus,  the  downtime  cost  function  can 
be  written 

fi(A-i.Hi)  -  FiKXiMi  *  (SU)iKXi  (24) 

where 

(1)  fi(Xi,M^)  is  the  cost  in  dollars  of  maintaining  subsystem  i 
with  repair  time  and  failure  rate  X^, 

(2)  Fi  is  the  cost  factor  representing  the  cost  per  hour  of  down¬ 
time  for  subsystem  1, 

(3)  (SU)j.  is  the  cost  factor  representing  the  fixed  cost  per  failure 
for  subsystem  i,  and 

(4)  K  is  the  total  system  operating  time,  in  hours,  over  the  life 
cycle  of  the  system. 

The  inventory  costs  are  associated  with  the  cost  of  procuring  spares, 
and  the  number  of  spares  needed.  For  example,  a  system  with  a  high  fail¬ 
ure  rate  would  require  many  more  spares  than  a  system  with  a  low  failure 
rate.  Although  for  very  low  failure  rate  items  the  sparing  cost  would 
not  be  a  linear  function  of  the  failure  rate,  for  a  large  number  of 
systems,  with  adequate  inventory  policies,  the  sparing  cost  would  gener¬ 
ally  be  a  linear  function  of  the  failure  rate.  Since  productioo  costs 
were  assumed  to  be  inversely  proportional  to  failure  rates  and  repair 
times,  it  would  be  reasonable  to  expect  that  the  oost  of  spares  would 
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also  be  inversely  proportional  to  repair  times  and  failure  rates.  Thus, 
the  sparing  cost  function  assumed  in  this  paper  is 


si(xi*Mi) 


SiKXi  _  SjK 

XiMx  "  Mi 


s1(Mi) 


(25) 


where 

(1)  sj(M^)  is  the  cost  in  dollars  of  sparing  subsystem  i  with 
repair  time  Mi, 

(2)  is  the  cost  factor  associated  with  sparing  subsystem  i, 
and 

(3)  K  is,  as  before,  the  total  system  operating  time,  in  hours, 
over  the  life  cycle. 

At  this  point,  the  problem  of  what  to  do  with  depot  maintenance 
costs  should  be  discussed.  First  it  should  be  noted  that  the  repair 
time,  as  discussed  to  this  point,  is  the  time  required  to  bring  the 
system  back  into  operating  condition  at  the  operational,  or  lower 
levels  of  maintenance.  It  may  or  may  not  have  anything  to  do  with 
depot  maintenance  time  which  could  be  termed  mean  time  to  overhaul  or 
rebuild.  If  the  cost  to  rebuild  a  unit  one  time  at  the  depot  level 
is  independent  of,  say,  remove  and  replace  time  and  reliability  at  the 
organizational  level,  then  depot  maintenance  costs  could  generally  be 
included  in  maintenance  set  up  costs.  However,  if  the  line  replaceable 
unit  is  made  larger  as  a  means  of  decreasing  repair  time,  depot  main¬ 
tenance  costs  would  depend  0:1  MTTR  as  previously  defined.  For  example, 
in  one  radio  set  the  line  replaceable  unit  may  be  an  IF  strip,  and  in 
the  second  case,  the  entire  receiver.  In  the  first  case  the  depot  would  , 
rebuild  the  IF  and  in  the  second  case,  the  entire  receiver.  In 


this  paper  it  is  assuaed  that  the  second  case  applies,  and  that  depot 
maintenance  costs  are  invorsely  proportional  to  both  repair  tine  and 
failure  rate.  Depot  maintenance  costs  are,  therefore,  included  in  the 
sparing  costs. 

Functions  describing  the  development  cost  of  increasing  the  avail¬ 
ability  from  a  level  below  the  requirement  to  the  required  level  should 
have  many  of  the  characteristics  of  the  production  cost  functions, 
except  that  if  the  optimal  allocation  Indicated  a  failure  rate  or  repair 
time  greater  than  that  already  achieved  no  saving  in  devalopment  cost 
would  be  realized  in  purposely  increasing  a  failure  rate  or  repair  time. 
Therefore,  the  development  cost  functions  chosen  for  this  paper  have  the 
form  of  the  curve  in  Figure  2. 


Xj_  allocated  failure  rate 
^  achieved  failure  rate 
allocated  repair  time 
achieved  repair  time 

Figure  2.  Development  Cost  Versus  Repair  Time  or  Failure  Rate 

This  curve  represents  &  characteristic  relationship  between  develop¬ 
ment  cost  and  repair  time  or  failure  rate  and  further,  it  indicates  that 


cost  is  not  increased  if  a  failure  rate  or  repair  tine  is  increased. 
Development  funds  already  spent  on  the  project  are  not  included  in  the 
nodel  since  the  purpose  of  the  allocation  is  to  find  the  optimal  set 
of  repair  tines  and  failure  rates  given  an  achieved  set  that  falls 
short  of  the  availability  goal. 

The  functions  used  to  describe  development  costa  for  improvement 
of  subsystem  i  are 


,  „  x  (DF)i  (DF)i  ,  „  p 

dfi(Xi)  -  ^ - jr-  <  *i* 

(26) 

dfi(Xi)  “Of  Xi  >  $1, 

(27) 

(mJi  (DM^ 

dml(Mi)  -  „  1  -  -3H1  »  «i  <  Mi. 

Mi 

(28) 

dmi(Mi)  -  0  »  Mi  >  I?i, 

(29) 

where 

(1)  dfi(A^)  is  the  cost  of  improving  the  achieved  failure  rate 

to 

(2)  d-ji(Mi)  is  the  cost  of  improving  the  achieved  repair  time  l?i 
to 

(3)  (DF)i  is  the  cost  factor  associated  with  the  difficulty  of 
decreasing  the  failure  rate  of  subsystem  i,  and 

(4)  (DM)i  is  the  cost  factor  associated  with  the  difficulty  of 
decreasing  the  repair  time  of  subsystem  i. 

The  dimensions  of  all  cost  factors  are  shown  below 
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PFp 

l/hr. 

1  PFi  =  $/hr. 

PMp 

hr. 

!  PMi  -  ($)(hr.) 

Fi(hr.) 

;  Fp  -  $/hr. 

SUj_hr. 

hr. 

;  SUi  =■  $ 

Sphr. 

hr. 

;  Si  =*  $ 

DFi 

1/hr. 

}  DFi  "  $/hr. 

DMi 

hr. 

>  DMi  -  ($)(hr.) 

Information  for  determining  the  hourly  maintenance  cost  must,  of 
course,  come  from  the  using  organization.  It  could  be  determined 
reasonably  from  p*st  data  and  should  include  training  costs  if  the 
technician  turnover  rate  is  high  enough  to  require  that  new  techni¬ 
cians  must  be  trained  throughout  the  system  life  cycle.  Maintenance 
set  up  costs  will  be  primarily  dependent  on  the  location  of  mainte¬ 
nance  technicians,  tools  and  test  equipment  relative  to  the  operating 
system,  and  the  cost  of  system  shut  off  and  start  up. 

As  implied  by  Equation  25,  the  cost  per  failure  is  given  by 


$/failure  =  ^ 

which  includes  the  cost  of  procuring  the  spare  and  carrying  it  in 
inventory.  Suppose  subsystem  i  has  the  following  characteristics, 

“  ,001/hr, 

“  1  hr. 


and  it  has  been  determined  that  the  average  inventory  cost  per  failure  i 


$/f all ure  -  $25 


then  Si  is  given  by 

Si  -  $25(.00l)  °  $.025. 

The  above  cost  factors,  hourly  maintenance,  maintenance  set  up 
and  sparing  costs,  should  all  be  weighted  by  an  expected  value  infla¬ 
tion  factor  and  a  present  worth  interest  factor.  They  may  also  be 
dependent  on  the  location  and/or  environment  of  the  operating  system. 

The  value  of  development  and  production  cost  factors  may  b°  depen¬ 
dent  on  the  level  of  availability  already  achieved  or  predicted  for  a 
given  design.  For  example,  a  proposed  system  design  may  have  a  pre¬ 
dicted  availability  of  .90  and  it  is  determined  that  a  .95  availability 
is  necessary.  In  this  case,  it  may  be  necessary  to  reevaluate  the 
development  and  production  cost  factors  since  different  technologies 
may  be  needed  for  such  an  improvement.  In  any  case,  the  total  devel¬ 
opment  cost  per  system  will  be  inversely  proportional  to  the  number  of 
systems  produced,  whereas  production  cost  per  system  will  be  relatively 
independent  of  the  number  of  systems  produced  if  production  set  up  costs 
are  included  in  development  costs  rather  than  production  costs. 

Since  the  mean  time  to  repair  is  typically  much  less  than  the  mean 
time  between  failures  the  failure  rate  development  cost  factor  should 
be  much  less  than  the  repair  time  development  cost  factor  if  develop¬ 
ment  efforts  required,  for  the  same  percentage  improvement,  are  similar. 
Thus,  given  a  subsystem  for  which  it  has  been  determined  that  to 
decrease  X.i  and  Mi  by  one  half  would  cost  $10,000  each,  and  the  parame¬ 
ters  of  the  subsystem  are  as  given  above,  that  is, 


Al  -  .OOl/hr 
Mj_  “  1  hr. 


then  the  cost  associated  with  each  would  be 

DF< 

•  Tooos 

$10,000  - 

o 

or, 

DFi  -  $5/hr. 

DMi  -  $5,000/hr. 

Production  cost  factors  would  be  determined  In  a  similar  manner, 
however  the  reason  for  production  costs  varying  with  repair  time  and 
failure  rate  would  be  due  to  such  factors  as  more  expensive  materials, 
more  elaborate  production  techniques,  or  tighter  quality  control* 

All  cost  factors  are  constants  that  are  determined  for  each  sub¬ 
system  based  on  existing  data,  past  experience,  familiarity  with  the 
equipment,  and  engineering  judgment.  The  importance  of  careful  esti¬ 
mates  for  these  factors  cannot  be  overemphasized. 

If  the  functions  previously  discussed  are  substituted  into  the 
total  cost  equation  obtained  in  Chapter  III,  the  following  equation  Is 
obtained i 


J3Fi  +  PFi  DMi  +  PMi  SiK 

f .  C .  **  U_— 


+  jjj~  ♦  SUiKXi 


DM. 


i  DF«  uii 

+  FjKXiMiJ  +  -  jr  +  1)  -  -  -jjr-  , 

Differentiating  with  respect  to  and  Q  yields 

dT^G.  „  TSDFi't'.,FFl)i.  SUlK  +  PlKMl  +  QMi  -  0 


(30) 


(31) 


i 
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for  i  **  1,  2,  . . . ,  n 


for  i  “  1,  2,  ... i  n  and 


^  *  F  +  1  -  0 


(32) 


(33) 


The  solution  to  this  set  of  equations  represents  the  set  of  improve¬ 
ments  in  subsystem  repair  times  and  failure  rates  which  minimize  the 
improvement  cost  while  achieving  the  availability  goal. 

Separating  the  variables  Ap  and  Mp  yields 

(SUiK)aXi  -  (FjK  +  Q)(DMi  +  PMp  +  SjKHf 

-  2  SUiK(DF1  +  PFi)X?  +  (DFi  +  PFp)3  ■*  0  (34) 

f*or  i  **  If  2  $  n 

(FiK  +  Q)a(DFi  +  PFi)Mi  -  (DM*  +  PMi  +  SiK)a(F1K  +  Q)M^ 

-  (DMi  +  PMi  +  SiK)aSUiK  ■=  0  (35) 

for  i  -  1,  2,  ...  n 

Although  general  solutions  to  fourth  order  polynomials  exist,  it  is 
unlikely  that  they  would  add  any  insight  to  the  solution  in  this  case, 

A  computer  program  is  listed  in  the  Appendix  which  was  developed  to 
solve  the  equations  by  the  Newton-Raphson  technique,  first  with  the 
constraint  inactive,  and  second,  with  the  constraint  active. 

If  the  hourly  maintenance  cost  for  each  subsystem  is  the  same, 


Equation  30  can  be  written  as 
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-DF<  +  PFi  DM1  +  PMi  +  SiK 
T.C.  -  +  -i - - - -  +  SUiKXj 

+  F  K  I  A^Mj_  +  Q(2-A^M^  —  -~r  i  1)  ■  ~jA  — 


where  F  is  equal  to  F^(  a  constant  for  all  1, 
The  above  equation  can  also  be  written 


pDFi  +  PFi  EKi  +  PMi  +  SiK  - 

r.c.  -  it- — jq- —  + - - +  syiKXiJ 


since,  when  the  constraint  is  active 


(36) 


2^1%  ■  jr  -  1« 

Here  ft  is  seen  that  under  the  conditions  of  identical  hourly  mainte¬ 
nance  costs  for  each  subsystem,  and  an  active  availability  constraint, 
the  optimal  values  of  Xi  and  Mi  are  independent  of  hourly  maintenance 
costs.  The  reasoning  is  that  an  active  availability  constraint  implies 
that  a  fixed  amount  of  maintenance  will  be  performed  and  if  all  hourly 
maintenance  costs  the  same  it  makes  no  difference  on  which  subsysteu 
it  la  performed, 

A  further  simplification  can  be  made  if  maintenance  set-up  costs 
are  negligible.  This  would  imply  that  the  maintenance  technicians  and 
necessary  repair  equipment  are  in  close  proximity  to  the  operating  sys¬ 
tem  and  also  that  shut-down  and  st-rt-up  costa  are  negligible.  If  this 
ia  the  case  Equation  36  can  be  written  as 

t.c.  -  g- °Fi  t 31  *  *£Uh 

•  Cl  m  J 


+ 


reZXiMi  +  Q(ZXiMi  -  JT  +  1) 


DFi  DMi 

Ti’Ti 


(37) 
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Differentiating  now  with  respect  to  A^,  M^,  and  Q  yields 


iliSa  +  (FK  +  QjMi  -  0 

3  XL  Aj  A 

for  1  *  lj  2|  a • i i  n 


iliSa 

3  Mi 


"(DMi 


+  PMi 

4 


SlK) 


+  (FK  +  Q)Ai  “  0 


(38) 


(39) 


for  i  -  1,  2,  n 

and 


~Jq~  -  jr  +  1  -  °. 

Solving  for  A^,  M^,  and  Q  yields 

,  ...  r  (»Fi  +  PFi)3 

h.  “  '“'(TIC”  T  QKWi  +  PMi  T  SiK)J  » 

u  .  r  (nil  +  PM1  I  si*)3-,l/3 

1  ‘"Tra  +  Q)(DF1  +  PF^ 

for  i  -  1,  2,  ....  n 


and 


-  1 


2/3 


,  ^  r  _ A  _ -f/J 

°  "  L  (DFi  +  PFi)(mi  +  PMi  +  S^K)“*  “  ^ 


(40) 


(42) 


Although  the  assumption  that  the  hourly  maintenance  cost  is  the  same  for 
each  subsystem  may  apply  to  a  large  variety  of  systems,  the  assumption 
that  maintenance  set-up  costs  are  zero  is  quite  restrictive.  Here 
again,  all  assumptions  must  be  carefully  investigated  before  choosing 
a  model. 

The  computer  program  listed  in  the  Appendix  will  solve  all  cases 
discussed  to  this  point.  If  hourly  maintenance  costs  are  constant  for 


all  subsystems  they  are  punched  the  sane  on  all  data  cards,  and  if 
maintenance  set  up  costs  are  zero  they  are  read  in  as  zero* 

If,  for  the  equations  given,  an  allooated  failure  rate  (X^)  or 
repair  time  (Mj,)  is  greater  than  that  already  achieved,  the  develop- 
aent  cost  factor,  DF^  or  DM^  respectively,  is  set  equal  to  zero  since 
a  zero  cost  is  associated  with  increasing  a  failure  rate  or  repair  tine. 
The  equations  are  then  resolved  for  the  sane  value  of  the  Lagrange 
multiplier.  If,  upon  resolving,  the  respective  failure  rate  or  repair 
tine  is  less  than  that  already  achieved,  the  optimal  value  is  identical 
to  the  value  already  achieved*  If.  however,  the  newly  allocated  failure 
rate  (repair  time)  is  greater  than  that  already  achieved,  the  optimal 
value  is  given  by  the  value  obtained  with  DFi(DMi)  -0,  This  can  be 
illustrated  an  follows i  the  development  cost  for  failure  rate  i  is 
given  by  Equations  26  and  27  which  are  repeated  here  for  convenience, 

.  DF<  DFi  . 

dfi(*i)  "  for  *1  <  4 

dfi(*i)  “  0  for  >  ^i  • 

The  first  equation  intercepts  the  zero  axis  at  X^  O  but  is  valid 
only  for  X^  <  ^,  thus,  the  development  cost  curve  is  continuous,  and 
the  total  cost  is  as  shown  in  Figure  3* 


TOTAL 

GOST 
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Curve  1,  Cost  curve  assuming  DFp  >  0  for  all  Xp 
Curve  2.  Cost  curve  assuming  DFi  =  0  for  all  Xp 
______  Actua'1  Cost  Curve 

Figure  3 «  Total  Cost  versus  Failure  Rate 

Suppose  £p  is  the  value  already  achieved.  Curve  1  applies  if  Xp  <  £p 
while  curve  2  applies  where  Xp  >  £p.  While  the  optimal  value  of  curve  1 
is  to  the  right  of  £p,  the  optimal  value  of  curve  2  lies  to  the  left  of 
£p.  However,  the  optimal  value  of  the  total  cost  curve  lies  at 
=  Xp. 

The  case  above  may  indicate  that  the  failure  rate  or  repair  time 
in  question  was  achieved  at  more  expense  than  necessary,  while  a 
cheaper,  more  advantageous  route  to  the  availability  requirement 
existed.  Or,  it  may  be  the  result  of  the  various  cost  factors  changing 
due  to  either  unforeseen  problems  or  technological  advances  during 
research  and  development. 

It  would  be  of  interest  to  obtain  the  optimal  failure  rates  and 
repair  times  for  a  system  that  has  no  specific  availability  requirement. 
Of  course,  in  this  case,  hourly  maintenance  costs  must  be  included, 
whether  or  not  It  is  the  same  for  each  subsystem,  since  there  is  no 
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constraint  on  the  amount  of  maintenance  to  be  performed.  The  computer 
program  in  the  Appendix  will  obtain  an  optimal  set  of  failure  rates  and 
repair  times  (the  first  set  in  the  print-out  under  the  heading  "The 
Lagrange  Multiplier  is  Zero")  however,  it  is  assumed  that  the  system 
will  operate  a  fixed  number  of  operating  hours. 

This  concludes  Chapter  IV.  The  solutions  obtained  in  this  chapter 
are  applied  to  an  example  problem  in  Chapter  V  and  the  computer  program 
developed  to  handle  the  calculations  is  discussed. 


CHAPTER  V 


APPLICATION  OF  ALLOCATION  TECHNIQUE 

An  example  problem  is  considered  in  this  chapter  for  a  system 
which  consists  of  only  five  subsystems,  although  the  model  and  the 
computer  program  can  be  used  to  handle  a  system  of  any  size. 

To  enable  solution  of  the  problem  by  the  program  listed  in  the 
Appendix  the  following  parameters  must  be  known i 

(1)  Development,  production,  sparing  and  maintenance 
cost  factors. 

(2)  Required  inherent  availability. 

(3)  Acceptable  availability  tolerance. 

(4)  Expected  total  operating  hours  of  thi  system. 

(5)  Predicted  values  of  failure  rates  and  repair  times 
currently  achieved. 

Determination  of  the  various  cost  factors  was  briefly  discussed 
in  Chapter  IV,  and  as  previously  stated, this  paper  is  not  concerned 
with  establishing  an  availability  requirement.  An  availability  toler¬ 
ance  is  necessary  to  ensure  that  the  computer  does  not  iterate  the 
solution  procedure  endlessly  if  it  cannot  converge  exactly  to  the 
required  availability  due  to  truncation  error.  For  example,  the  avail¬ 
ability  requirement  may  be  stated  as 

A  ■  ,95  &  .001 
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Determination  of  the  expected  total  operating  nours  of  the  system 
may  be  quite  difficult  since  it  is  dependent  on  administrative  and 
logistics  downtime  which  may  vary  with  the  using  organization  and  the 
location  of  the  system.  Total  operating  time  will  also  be  dependent 
on  the  amount  of  preventive  maintenance  performed. 

To  illustrate  how  this  figure  could  be  determined,  suppose  that 
the  administrative  and  logistics  dslay  time  is,  typically,  some  per¬ 
centage  of  the  hourly  maintenance  time,  and  that  the  system  is  either 
operating  or  it  is  down  for  maintenance.  The  total  calendar  time  can 
then  be  expressed  by  the  following  equation* 

T  -.TAq  +[(1/Aj;  -  l)P  +  (l/Aj  -  1^  +  *  TAq 

then, 

+ "(i/Ai~lTO'+"p) 

,  where 

Aq  is  the  percentage  of  the  total  time  that  the  system  is  operating, 

(l/Aj  -  1)  represents  the  amount  of  corrective  maintenance  per  hour 
of  operating  time, 

P  is  the  amount  of  administrative  and  logistics  delay  time  per 
hour  of  corrective  maintenance, 

Aj  is  the  inherent  availability  as  previously  defined, 

is  the  amount  of  preventive  maintenance  required  per  hour  of 
operating  time,  and 

T  is  the  total  calendar  time  in  hours. 

Thus,  the  total  operating  time  of  the  system  is  given  by 


Operating  time  "  TA<» 


It  was  assumed  that  no  administrative  or  logistics  delay  time  is 
associated  with  preventive  maintenance!  however,  tb-:  total  time 
aquation  could  easily  he  written  to  include  such  a  term* 
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If  an  initial  allocation  is  being  performed,  that  is,  no  values 


for  achieved  failure  rates  and  repair  times  exist,  the  last  two  terms 
DFi  EMi  9 

of  Equation  30,  -jp  and  *-j-,  are  zero,  since  and  %  axe  virtually 

infinite  and  some  nonzero  development  t  will  be  associated  with 

obtaining  finite  repair  tlmeB  and  failure  rates.  When  using  the  program 

listed  in  the  Appendix  it  is  necessary  to  enter  unrealistically  high 

values  for  the  achieved  failure  rates  and  repair  times  when  performing 

DF^  LMi 

an  initial  allocation.  When  this  is  done  the  terms  and  —5-  will  be 

Mi 

Y 

very  small  and  can  be  neglected  in  the  total  cost  equation. 

The  input  format  and  sequence  of  all  cost  factors,  and  the  loca¬ 
tion  of  the  parameters  discussed  above  is  given  In  the  Appendix. 

The  example  problem  solved  is  for  the  case  in  which  a  system  has 
been  proposed  that  falls  short  of  the  availability  goal.  Predictions 
for  the  failure  rates  and  repair  times  are  listed  in  Table  1.  The 
repair  time  includes  all  active  portions  of  downtime.  For  example, 
these  might  consist  of  fault  isolation,  removal,  repair  and/or 
replacement,  and  check  out. 
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Table  1.  Subsystem  Repair  Times  and  Failure  Rates 


Subsystem 

Repair  Time(hr.) 

Failure  Rate(failure/hr.) 

1 

15 

.0019 

2 

8 

.0051 

3 

10 

.0034 

4 

10 

.0034 

5 

10 

.0034 

Availability  -  .8538 

The  cost  factors  assumed  for  this  example  are  listed  in  Table  2. 


Table  2.  Subsystem  Cost  Factors 


Subsystem 

DM($-hr.) 

PM($-hr.) 

DF($lhr.) 

PF($lhr.) 

1 

17,000  ’ 

17,000 

25.00 

25.00 

2 

12,000 

12,000 

18.00 

*3,00 

3 

20,000 

20,000 

10.00 

10.00 

4 

20,000 

20,000 

5.00 

5.00 

5 

30,000 

30,000 

10.00 

10.00 

S($) 

su($) 

Fi($lhr.) 

1 

.10 

2.00 

8.00 

2 

.06 

2.00 

6.00 

3 

.05 

2.00 

8.00 

4 

.08 

2.00 

10.00 

5 

.12 

2.00 

12.00 

It  is  desired  to  achieve  an  availability  of  .95  at  the  least  possible 
cost. 

The  optimal  allocations  for  this  system,  found  by  Messer  [7]  from 
considering  only  development  costs,  are  listed  in  Table  3* 


Table  3.  Subsystem  Allocations  Considering  Development  Costs  Only 


Subsystem 

Repair  Time(hr.) 

Failure  Rate (f allures/ hr.) 

1 

3.10 

.o0190 

2 

2.?4 

.00411 

3 

4.69 

.00234 

4 

5.90 

.00146 

5 

A  < r tm  4  1  V4  1  4  — 

6.14 

OCOQ 

.00205 

Availability  -  .9529 


Comparing  these  results  to  the  repair  times  and  failure  rates  in  Table  1, 
it  is  seen  that  the  first  failure  rate  has  remained  the  same  while  all 
other  parameters  have  decreased. 

Substituting  all  of  the  cost  factors  listed  in  Table  2  into 
Equations  33.  34,  and  35  yields  the  results  listed  in  Table  4-, 


Table  4,  Subsystem  Allocations  Considering  All  Cost  Factors 
Subsystem  Repair  Time(hr.)  Failure  Rate(failures/hr. ) 


1 

3.9955 

.002741 

2 

2.8416 

.003988 

3 

4.6910 

.002265 

4 

5.8693 

.001400 

5 

6.0118 

.001913 

Availability  -  .9500 

Comparing  Table  4  and  Table  3  it  is  seen  that  the  first  failure  rate 
should  deliberately  be  increased.  Messer  used  the  same  development 
cost  function  that  is  used  here,  namely  that  the  development  cost  for 
increasing  a  failure  rate  or  repair  time  is  zero,  so  that,  if  develop¬ 
ment  costs  only  are  considered,  no  repair  time  or  failure  rate  could 
increase  in  the  optimal  solution.  However,  when  all  cost  factors  are 
considered  it  may  be  desirable  to  increase  one  or  more  failure  rates  or 
repair  times.  If  the  allocation  procedure  indicates  that  a  failure  rate 
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or  repair  tine  should  be  increased,  no  savings  In  development  costs  can 
be  realized  since  development  funds  to  achieve  the  present  levels  of 
repair  tines  and  failure  rates  have  already  been  spent.  However,  a 
savings  in  production  costs  will  be  realized  for  the  case  where  the 
solution  indicates  an  increasing  failure  rate,  and  a  savings  in  produc¬ 
tion  and  sparing  costs  can  be  realized  for  the  case  of  an  increasing 
repair  time. 

The  results  of  the  allocation  technique  are  guidelines  to  be  taken 
into  consideration,  not  hard  and  fast  rules.  Also,  it  will  never  be 
more  accurate  than  the  estimates  of  the  cost  factors  involved,  no’-  '-ore 
valid  than  the  assumptions  made  in  deriving  the  model. 

The  total  additional  cost  is  also  evaluated  at  each  iteration  of 
the  program,  however,  it  must  be  cautioned,  that  one  should  not  compare 
the  costs  listed  as  a  means  of  determining  the  "price"  of  additional 
availability.  As  stated  before,  the  allocations  are  made  under  the 
assumption  that  the  system  will  operate  a  fixed  number  of  hours j  there- 

s' 

*fore,  if  the  availabilities  listed  are  significantly  different  from 
the  required  availability,  system  lifetime  must  also  be  significantly 
different  from  the  planned  lifetime  of  the  system. 

Some  estimates  of  the  sensitivity  of  failure  rates  and  repair  times 
to  the  various  cost  factors  are  obtained  in  Chapter  VI. 


CHAPTER  VI 


SENSITIVITY  ESTIMATION 

Knowledge  of  the  sensitivity  of  failure  rates  and  repair  times  to 
the  various  parameters  discussed  earlier  is  important  when  determining 
how  much  effort  should  be  devoted  to  estimating  each  parameter.  If 
repair  times  and  failure  rates  are  relatively  insensitive  to  one  param¬ 
eter  as  opposed  to  another  it  is  obvious  that  more  effort  should  be 
alloted  to  obtaining  an  accurate  estimate  of  the  more  sensitive 
parameter. 

In  this  chapter  an  attempt  is  made  to  obtain  estimates  of  failure 
rate  and  repair  time  sensitivity  by  varying  certain  parameters  and 
observing  the  effects  of  their  variation.  This  example  consists  of  a 
system  with  five  subsystems  having  the  same  cost  factors.  They  are  as 
shown  in  Table  5» 

The  number  of  operating  hours  Is 

K  =  20,000, 

and  the  availability  constraint  is 

Aj  -  .95  ±  .0001. 

The  results  of  this  allocation  will,  of  course,  have  all  Xi's  equal  and 
all  M^'s  equal.  They  are 

Xi  -  X2  “  ...  -  X5  -  .00285368 

Mi  -  M2  -  ...  -  M5  -  3.o872. 

U0 


41 


Table  5 

Cost  Factors  for  Sensitivity  Estimation 


Subsystem 

DFi 

PFi 

DMi 

PMi 

1 

10 

10 

10,000 

10,000 

2 

10 

10 

10,000 

10,000 

3 

10 

10 

10,000 

10,000 

4 

10 

10  ■ 

10,000 

10,000 

5 

10 

10 

10,000 

10,000 

Subsystem 

Si 

SUi 

Fi 

0 

0 

Mi 

1 

.25 

4.0 

12.0 

99999. 

99999999 

2 

.25 

4.0 

12.0 

99999. 

99999999 

3 

.25 

4.0 

12.0 

99999. 

99999999 

4 

.25 

4.0 

12.0 

99999. 

99999999 

5 

.25 

4.0 

12.0 

99999. 

99999999 

The  availability  is 

Ai  -  .9500, 


and  the  total  cost  is 


T.C.  =  $82,711.25  per  system. 


By  recalling  Equation  30, 


DFi  +  PFi  DMi  +  PMi  +  SiK 

T.C.  =  E[ - r -  + - - -  +  SUjKXi  +  PiKAiM] 


it  is  obvious  that  the  sensitivity  of  repair  times  and  failure  rates  will 
not  be  dependent  on  a  percentage  change  in,  say,  DFi  alone,  but.  on  the 
percentage  change  of  the  sum  DFi  +  PFi.  Likewise,  it  will  be  dependent 
on  the  percentage  change  of  the  sum  DMi  +  PMi  +  SiK.  In  the  estimation;-, 
that  follow  only  one  parameter  at  a  time  was  varied  while  all  others  were 
as  listed  in  Table  5«  All  variations  were  made  on  subsystem  5*  The 
results  of  these  estimates  are  shown  in  Table  6. 
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Table  6 


Sensitivity  to  Cost  Factors 


Xl  -  14 


DF5-H>F^-40  .002785 

DF5+PF5-IO  .002912 

DMc-H?M5+S5K-50,000  .002783 

mi4fHc*Sdi-lZ,500  .002913 

SU5-8  .002856 

SU5-2  .002852 

F5-24  .002912 

F5-6  .992817 


b 

Mj  -  M4 

m5 

.004398 

3.5933 

2.8596 

.001843 

3.7678 

4.7361 

.002226 

3.5902 

5.6773 

.003624 

3.7699 

2.3903 

.992797 

3.6907 

3.7296 

.002883 

3.6852 

3.6652 

.002617 

3.7675 

3.3628 

.003000 

3.6369 

3.8896 

DF  5+PF  5“40 
DF5+PF5-IO 
DM5+PM5+S5K-50 ,000 
DM5+PM5+S5K-12.000 
SU5-8 

SU5-2 

F5-24 

F5-6 


Cost  2 

Change  in 
Parameter 

%  Change 
in  Cost  1 

83,710.86 

doubled 

6.7 

83.853.56 

halved 

-5.2 

83,681.38 

doubled 

6.5 

83,791.15 

halved 

-5.0 

82,729.18 

doubled 

.28 

82,711.80 

halved 

-.14 

82,879.25 

doubled 

2.8 

82,748.78 

halved 

-1.6 

Cost  1 

88,257.56 

78,424.88 

88,084.31 

78,578.19 

82,942.69 

82,595.50 

84,990.75 

81,365.06 


%  Change 
in  Cost  2 

1.2 

1.4 

1.2 

1.3 

.022 

.00067 

.20 

.045 


Cost  1  represents  the  true  cost  per  system  if  the  cost  factor  listed  In 
the  left  hand  column  actually  has  the  value  stated,  whereas  Cost  2 
represents  the  true  cost  per  system  if  the  cost  factor  in  the  left  hand 
column  is  estimated  to  have  the  value  given  and  has  its  actual  value  as 
listed  in  Table  5« 

Although  no  general  conclusions  can  be  drawn  from  Table  6,  since 
specific  numbers  were  assumed,  the  table  indicates  that  the  optimal 
repair  times  and  failure  rates,  as  well  as  the  optimal  cost,  are  rela¬ 
tively  insensitive  to  the  maintenance  cost  factors,  SUp  and  Fi.  This 
should  be  expected  since,  when  an  availability  constraint  exists,  the 
total  maintenance  performed  will  be  forced  to  be  relatively  small  com¬ 
pared  to  the  optimal  unconstrained  value,  unless,  of  course,  the  opti¬ 
mal  value  is  nearly  equal  to  the  constrained  value. 

The  sensitivity  with  respect  to  the  number  of  operating  hours  was 
obtained  using  the  cost  factors  of  Table  5»  and  varying  K,  The  results 
are  listed  in  Table  7. 


Table  7 


Sensitivity  to  the  Number  of  Operating  Hours 


*i 

Mi 

Cost  per  System 

20,000 

.002854 

3.6872 

$82,711.25 

10,000 

.003032 

3.4727 

$72,316.50 

40,000 

.002577 

4.0830 

$102,854.60 

‘  I  Table  7  indicates  that  as  the  number  of  operating  hours  increase 


th:.  optimal  failure  rates  decrease  while  the  optimal  repair  times 
in.rease.  The  reasons  for  this  are  that  the  total  sparing  cost  is 
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proportional  to  K  and  inversely  proportional  ro  and  the  total  main¬ 
tenance  set  up  cost  is  proportional  to  both  K  and  A^. 

Finally,  the  sensitivity  with  respect  to  the  number  of  systems 
produced  was  estimated  by  varying  the  development  cost  factors. 

Reducing  DFi  and  DMi  by  one  half  implies  that  twice  as  many  systems 
as  represented  by  Table  5  will  be  produced  while  increasing  DF^  and 
DMi  by  a  factor  of  two  implies  that  half  as  many  systems  will  be  pro 
duced.  All  other  cost  factors  remained  the  same.  The  results  of  these 
allocations  are  shown  in  Table  8, 

Table  8 


Sensitivity  to  Number  of 

Systems  Produced 

*i 

Mi 

Cost  per  System 

DFi 

DM-. 

*=  10 
*  10,000 

.00285368 

3.6872 

$  82,711.25 

DFi 

DMi 

Q  Q 

V_n 

O 

O 

O 

.00275190 

3.8236 

$  67,13^.^ 

DFi 

DMi 

=*  20 

-  20,000 

.00296773 

3.5^56 

$113,713.80 

In  all  cases  the  availability  is 

Aj  -  .9500 

Table  8  shows  that  the  optimal  repair  times  and  failure  rates  are 
also  relatively  insensitive  to  the  number  cf  systems  produced j  however, 
the  availability  constraint  is  also  active  here  and  the  optimal  alloca¬ 
tions  would  not  be  expected  to  change  significantly  when  all  development 
cost  factors  change  by  the  same  percentage.  However,  as  stated  before, 
the  change  in  or  Mi  is  not  sensitive  to  DFi  or  DMi  per  se,  but  to 


the  sums  DF^  +  PFj.  and  DM^  +  PM^  +  S^K.  Thus,  the  sajne  percentage 
change  in  two  development  cost  factors  may  cause  two  very  different 
changes  in  the  sums  stated  depending  on  the  relative  magnitude  of  the 
development  cost  factors  to  the  other  cost  factors.  Under  these  condi¬ 
tions  the  optimal  allocations  would  not  be  expected  to  be  insensitive 
to  the  number  of  systems  produced. 

In  the  following  chapter,  the  results  of  this  paper  are  summarized 
and  recommendations  for  further  study  are  given. 


CHAPTER  VII 


SUMMARY  AND  CONCLUSIONS 

A  method  for  the  cost-based  allocation  of  repair  times  and  failure 
rates  has  been  developed.  The  solution  is  based  on  the  minimization  of 
life  cycle  cost  subject  to  the  constraint  of  an  inherent  availability 
requirement, 

"is  allocation  technique  is  applicable  to  systems  in  which  all  sub¬ 
systems  exhibit  a  constant  failure  rate  and  failures  occur  independently. 
The  removal  of  these  assumptions  would  generalize  the  allocation  proce¬ 
dure  and  certainly  make  it  more  realistic.  However,  with  consideration 
of  tine  dependent  failure  rates,  analytic  solution  techniques  are  usually 
infeasible,  if  possible  at  all. 

It  is  also  assumed  that  the  system  can  be  described  by  a  series 
model;  that  is,  all  subsystems  are  necessary  for  proper  system  operation, 
and  that  when  a  failure  occurs  in  one  subsystem,  all  subsystems  are 
turned  off  for  the  duration  of  maintenance  downtime.  Extensions  to 
other  models  have  not  been  considered  in  this  paper,  although  it  would 
appear  feasible  to  solve  the  allocation  problem  for  cases  in  which 
definite  system  models  could  be  developed.  The  only  change  necessary 
is  that  of  formulating  the  availability  constraint. 

The  usefulness  of  this  technique  lies  in  its  ability  to  determine 
the  set  of  repair  times  and  failure  rate  improvements  that  result  in  the 
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minimum  life  cycle  cost  to  achieve  a  specified  availability.  Another 
useful  feature  of  the  technique,  as  applied  in  the  computer  program,  is 
that  the  optimal  repair  time  and  failure  rate  allocations  are  obtained 


for  the  case  where  the  availability  constraint  is  inactive.  However, 
a  different  system  lifetime  is  implied  if  the  resulting  availability  is 
significantly  different  than  that  of  the  constraint. 

It  was  determined  that  if  hourly  maintenance  costs  are  identical 
for  each  subsystem  then  the  optimal  repair  time  and  lure  rate  alloca¬ 
tions  are  independent  of  hourly  maintenance  costs  when  the  availability 
constraint  is  active.  It  should  be  noted,  however,  that  this  was  made 
possible  by  the  way  in  which  the  availability  constraint  was  formulated. 
No  provision  has  been  made  in  this  paper  to  account  for  the  effects 
of  various  spares  provisioning  policies,  manpower  levels  and  repair 
facilities  on  the  system  and  subsystem  repair  times.  Hell -designed 
systems  can  have  disastrous  field  results  when  accompanied  by  ill- 
planned  logistics  and  administration.  Therefore,  it  is  necessary  to  con¬ 
sider  the  effects  of  these  policies  on  total  Bystem  operation.  Even 
f  after  a  basic  allocation  urogram  has  been  instigated,  these  feu: tors 
should  be  taken  into  consideration  before  indicated  changes  are  made 
in  the  system. 

Specific  work  could  also  be  accomplished  for  cases  where  a  maximum 
allowable  repair  time  exists.  It  is  possible  that  a  system  could  meet 
the  availability  requirement  and  still  experience  a  failure  that  would 
result  in  a  long  downtime  that  could  produce  disastrous  results.  For 
this  reason  it  would  be  desirable  to  have  an  allocation  technique  which 
would  ensure  that  maximum  allowable  repair  times  were  not  exceeded.  For 


this  problem,  it  is  necessary  to  consider  repair  time  distributions  so 
that  variations  about  the  mean  can  be  analyzed. 

It  would  also  be  desirable,  in  eases  where  high  availability  is  not 
essential,  to  allocate  repair  times  and  failure  rates  subject  to  a  total 
production  constraint j  that  is,  allocate  the  parameters  so  that  a  mini¬ 
mum  required  amount  of  production  could  be  achieved  over  a  fixed  period 
of  time  at  a  minimum  cost.  Essential  to  this  problem  is  formulating  a 
production  constraint  for  a  given  system  configuration. 

In  conclusion,  the  results  of  this  paper  should  be  of  interest  to 
equipment  designers  and  users.  Vfhere  it  is  necessary  to  achieve  a 
required  inherent  availability,  this  allocation  technique  will  direct 
attention  to  the  specific  parameters  that  need  to  be  improved.  Also, 
forcing  the  contractor  to  consider  all  cost  factors  involved  will  pro¬ 
vide  the  user  with  valuable  information  on  logistics  and  administrative 
planning  even  if  some  r ?  the  cost  Information  must  cone  from  the  user. 


APPENDIX 


COMPUTER  PROGRAM  FOR  THE  ALLOCATION  PROBLEM 

The  purpose  of  the  program  listed  here  is  to  provide  a  computer¬ 
ized  technique  for  solving  the  availability  allocation  problem.  The 
program  was  run  on  the  IBM  360-651  however,  it  could  be  run  on  the 
IBM  1130  if  the  statements  in  lines  80  and  81  are  replaced  by  the 
following  statements! 

XF(Q)  91,  90,  91 

90  IF(A-AA)  91,  55,  55 

91  if(abs(a-aa)~acc)  55,  55,  1? 

The  program  requires  less  than  8,000  bytes  of  memory  for  a  system  com¬ 
posed  of  five  subsystems.  To  rewrite  the  program  for  a  larger  system 
it  is  necessary  to  make  two  changes  in  the  program.  These  are 

(1)  The  variables  in  the  DIMENSION  statement,  line  2,  must 

be  dimensioned  as  large  as  the  number  of  subsystems  in  the 
system,  and 

(2)  The  number  of  subsystems,  N,  must  be  stated  in  line  5» 

The  significant  variables  of  the  program  are  listed  below, 

followed  by  a  program  flowchart. 

Variables 

DF(l)  Failure  rate  development  cost  factor  of  ith  subsystem 
DM(l)  Repair  time  development  cost  factor  of  ith  subsystem 
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PF(l)  Failure  rate  production  cost  factor  of  ith  subsystem 

PM(l)  Repair  time  production  cost  factor  of  ith  subsystem 

S(l)  Sparing  cost  factor  of  ith  subsystem 

SU(l)  Maintenance  set-up  cost  factor  for  subsystem  i 

F(l)  Hourly  maintenance  cost  fate  tor  for  subsystem  i 

XF(l)  Achieved  failure  rate  of  subsystem  i 

XM(l)  Achieved  repair  time  of  subsystem  i 

AA  Required  availability 

A  Allocated  availability 

ACC  Availability  tolerance 

N  Number  of  subsystems 

K  Total  number  of  system  operating  hours 

Q  Lagrange  multiplier 

XFN'(l)  Allocated  failure  rate  of  subsystem  i 
XKN(l)  Allocated  repair  time  of  subsystem  i 
AU,  A3,  A2,  Al,  AO  Coefficients  of  failure  rate  polynomial 

B^,  33,  B2,  Bl,  BO  Coefficients  of  repair  time  polynomial 

JF(l)  Coefficient  of  failure  rate  development  cost  factor 
(value  is  0  or  1) 

JM(l)  Coefficient  of  repair  time  development  cost  factor 
(value  is  0  or  1) 

If  the  Newton  Raphson  technique  does  not  converge  an  error  message 
is  printed  out  prior  to  the  failure  rate  and  repair  time  allocations 
and  the  allocations  will  be  incorrect.  If  convergence  was  obtained  no 
error  message  will  be  printed  out.  The  first  page  of  the  printout  con¬ 
tains  the  cost  factors  in  the  following  orderi 
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DFi  PFi  DMi  PMi  Si  SUi 
All  other  pages  of  the  printout  contain,  in  the  following  order, 
The  value  cf  the  Lagrange  multiplier. 

The  failure  rate  repair  time  allocations, 

The  inherent  availability, 

The  total  cost  in  dollars,  and 
A  message  that  the  cost  is  optimal  for  the  given 
availability  if  all  roots  converged. 


The  data  format  is  as  shown  on  pages  53  and  54.  The  data  input 
sequence  is  CARD  1^,  CARD  2^ ,  CARD  I2,  CARD  2^,  *##,  CARD  1^,  CARD  2^, 


INPUT  FORMAT  CARD  1* 


INPUT  FORMAT  CARD 
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Presently  achieved  repair  time 
F  10.8 


Presently  achieved  failure  rate 
F  10.8 


Hourly  maintenance  cost  factor 
F  10.2 


Sparing  Cost  factor 
F  10.2 


Repair  time  production  cost  factor 
F  10.2 


Repair  time  development  cost  factor 


F  10.2 


Failure  rate  production  cost  factor 
F  10.2 


Failure  rate  development  cost  factor 
F  10.2 


MAIN  PROGRAM  FLOWCHART  (CON'T.) 


SUBROUTINE  FLOWCHART 


NilCN 


int  .*  1u  Sit  II-  Ini  OPflMAL  SULUl  IbN  lHAI  ALHfcAGV  ACMltVfcU 
UG  Jl  I *t  .N 


FUNCTION  FuMCN.C  >'Ci,  Ci.CO.il 


^GQUO.UO  2UCCU.C0  U« CdOOG  i.OCOOC  lO.OOOOU 


f  A  l  L  o  M  t  «  A  f  t 


MUUMt  »*ft  Hi  P  A  l  ri  Tint 


fAiwu*t  Hilt  AfcPAJM  1 l*C 
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